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Rounding of Phase Transitions in Cylindrical Pores 
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Phase transitions of systems confined in long cylindrical pores (capillary condensation, wetting, 
crystallization, etc.) are intrinsically not sharply defined but rounded. The finite size of the cross 
section causes destruction of long range order along the pore axis by spontaneous nucleation of do- 
main walls. This rounding is analyzed for two models (Ising/lattice gas and Asakura-Oosawa model 
for colloid-polymer mixtures) by Monte Carlo simulations and interpreted by a phenomenological 
theory. We show that characteristic differences between the behavior of pores of finite length and 
infinitely long pores occur. In pores of finite length a rounded transition occurs first, from phase 
coexistence between two states towards a multi-domain configuration. A second transition to the 
axially homogeneous phase follows near pore criticality. 

PACS numbers: 64.75Jk, 64.60.an, 05.70Fh, 02.70Tt 
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Fluids and fluid mixtures in nano- and microporous 
materials (pore diameters from 1 nm to 150 nm) play 
important roles in various industries (extracting oil and 
gas from porous rocks; use as catalysts or for mixture 
separation in the chemical and pharmaceutical industry; 
nanofluidic devices, etc.) [l|-l3j]. The interplay of finite 
size and surface effects strongly modifies the phase be- 
havior of such confined fluids [l|, |3|— Il9l] in comparison with 
the bulk. The vapor to liquid transition is shifted ( "capil- 
lary condensation" has well as critical points [3j, U, |9|, [l2| . 
Effects of wetting [20j on phase coexistence give rise to 
interesting patterns (plugs versus capsules versus tube 
structures etc. [7|). However, although various phase 
diagrams (different from the bulk) have been proposed 
(e.g. [l|, y, III, ls|, [Pa LUa LLZ|), many aspects hitherto arc 
not well understood. E.g., the "critical point" where ad- 
sorption/ desorption hysteresis vanishes seems to be sys- 
tematically lower than the critical temperature where the 
density difference between the vapor-like and liquid-like 
states vanishes [12J, in contrast to what theories have 
predicted [14j. 

However, a crucial aspect (stressed only in a few pio- 
neering studies I3|, p , and in the context of Ising/lattice 
gas models [2ll423| ) is the rounding of all transitions, 
caused by the quasi-one-dimensional character of a fluid 
in a long cylindrical pore with cross-sectional radius R. 
With the current progress of producing pores of well- 
controlled diameter varying from the nanoscale (carbon 
nanotubes 231425 1) to arrays of pores in silicon wafers 
[2q . up to 150 nm wide and of well-controlled length, 
experiments become feasible which are not plagued by 
effects of random disorder, which occur in porous glasses 
[l|, |2g. Thus, it is important to understand the phase 
transitions in pores more precisely, considering both the 
radius R and the length L of the pore as variables (the 
important role of L has so far been largely disregarded) . 
In the present Letter, we elucidate the rounding of vapor- 
liquid type transitions in cylindrical pores, based on 
Monte Carlo simulations of two generic models and a 



phenomenological theory. We show that, even in the ab- 
sence of precursors of wetting, two rounded transitions 
occur. Near the pore critical temperature at the pressure 
where vapor and liquid in the pore coexist, a rounded 
transition occurs from a axially homogenous state to a 
multi-domain configuration, where vapor-like and liquid- 
like domains alternate. The properties in this region de- 
pend strongly on R but not at all on L. In contrast, at 
lower temperatures the system makes a transition, where 
the full capillary is either in a vapor-like or a liquid- 
like state. The location of this transition depends on 
L, and the vapor to liquid transition is accompanied by 
a pronounced hysteresis. We also show that the effective 
(size-dependent) free energy exhibits well-defined spin- 
odals (as a finite size effect), but they do not control 
dynamics. Nucleation of domain walls becomes domi- 
nant when their free energy cost are small (of order of a 
few fcgT, T being the temperature; henceforth kg = 1). 
This domain wall nucleation explains why the hysteresis 
disappears far below the capillary critical region for small 
pores. 

The simplest model that already shows some of these 
effects is the two-dimensional (2-d) Ising model on the 
square lattice in the geometry of D x L strips with peri- 
odic boundary conditions in both directions 2lM23l. |28| . 
While this 2-d model lacks the surface effects due to 
the walls of real 3-d pores, it exhibits already the dis- 
appearance of hysteresis far below pore criticality, since 
the condition L >> D suffices to stabilize multi-domain 
states, as we will show below. Spins Si = ±1 at lat- 
tice sites i interact with their nearest neighbors with an 
energy J = 1, and an external field H. We apply the 
standard single spin-flip Monte Carlo algorithm [29( and 
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record the magnetization M = ^ Si/N (N = LD; the 

i=l 

lattice spacing being the unit of length) as a function 
of H at various T. We start out with all spins up and 
H = 0.05. The system runs for a "time" t = 2 • 10 6 
Monte Carlo steps per spin (MCS). Then we decrease H 



in steps of AH = 0.001, and run the simulation at each 
field for the same time, until we reach H = —0.05. Af- 
terwards, we reverse the process and increase the field 
stepwise by AH until we are back at H = 0.05. The 
width of the resulting hysteresis loops (Fig. [Th) strongly 
decreases with increasing T and for T > T n (L 7 D), the 
"hysteresis critical point", a hysteresis is no longer ob- 
served. However, when we record the probability dis- 
tribution P(M) for H = with the Wolff cluster algo- 
rithm [29j, |30j we observe that P(M) still exhibits peaks 
very close to the (exactly known |31| ) spontaneous mag- 
netization M s at temperatures T > Tq(L,D). While for 
L = 480, D = 10 these peaks can be followed up to about 
T = 2.1, the peaks occur up to about the critical temper- 
ature for D = 10 and smaller L, e.g. L = 40. However, 
at Tq(L,D) an important change also occurs in P(M): 
while for T < T Q {L,D) for a wide range of M P{M) is 
strictly independent of M (corresponding to a slab con- 
figuration which contains exactly two non-interacting in- 
terfaces [32j), for T > Tq(L,D) a third broad peak in 
P(M) appears at M = 0. An examination of snapshot 
pictures of the system (Fig. [Ud) reveals that this 3 rd peak 
is due to multi-domain configurations 28ll33ll34|. 

Such multi-domain configurations can in fact be pre- 
dicted when one computes the correlation length £ along 
the strip (Fig. UK) by transfer matrix methods (for L — > 
oo 0) or Monte Carlo (for L > f [H |H|). The latter 
estimates were extracted from the wave vector-dependent 
susceptibility x(fc) for k = k m i B = (2ir/L, 0) 
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and agree perfectly with the exact results. Thus, for 
very long strips the statistical errors arc also well under 
control. 

This correlation length below criticality (where well- 
defined domains exist) just measures the typical dis- 
tance between domain walls along the strip. The ap- 
proximation based on the (exactly known |35j) interfa- 
cial free energy a, £ w cxp(Da/T) becomes only accu- 
rate when £ > 10 5 , i.e. at temperatures much lower 
than those of interest for Fig. Jl This simply repre- 
sents the well-known argument |36| that long-range or- 
der in quasi-one-dimcnsional systems is destroyed due to 
the entropy gain of putting interfaces into the system. 
The free energy difference (relative to the single-domain 
state) for a state with n (non-interacting) interfaces is 
F = nF{ n t + nT\n(n/eL), where the total free energy 
cost of one interface is given by Fi nt = Da. 

The occurrence of the central peak at T near Tq(L, D) 
means that when T is raised at H = there is a transition 
from nonzero (\M\) for T < Tq(L, D) to a state with no 
order ((|M|) < M s ) for T > T a (L,D). Wc characterize 
this transition by the weight of the central peak of P(M), 
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defined as W = J P(M)dM/ j P(M)dM where ±m 
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FIG. 1: (a) Magnetization of Ising strips for L = 480, D = 
10 plotted vs. field H at T = 1.5B, 1.6A, 1.7«, 1.8*, 1.9 and 
2.0. Runs with decreasing H are shown as full curves, with 
increasing H as broken curves. A detailed analysis shows 
that the hysteresis disappears at To = 1.9 in this case, (b) 
Distribution P(M) vs. M for H = and T = 1.8 - 2.2 from 
bottom to top at M = 0. The inset is a typical snapshot at 
T — 2.1 containing multiple domains stretched in y-direction 
by a factor ~ 4. 



are the locations of the minima of P(M). Fig. [5] shows 
that the "equal weight" rule (first order transitions from 
one state to another state occur when the weights of the 
two states are W = 1/2) roughly corresponds to the con- 
dition £ ss L/3. With increasing L the transition gets 
shifted to lower temperature and also gets sharper. Since 
W w 0.1 for £ = L and W « 0.9 for £ = L/9, we use 
£ sa exp(£>cr/T) w exp(2D/T) for low T for L -> oo to 
estimate both the location of the transition and its width 
AT. 



T (L,D) 



2D 



ln(L/3)' 



L 



AT 



In 3 



oo, 



T (L,D) ln(L/ 



Finally, defining a barrier AF from P{M) as AF = 
Tln[P(M max )/P(m)] we see (Fig. [3ji) why the transi- 
tion at Tq(L, D) is related to the vanishing of hysteresis: 
AF(T) — > at a temperature where M max and m merge, 
which is close to To + AT where W — )• 1. Actually, the 
hysteresis already vanishes when AF/T w 10 since then 




FIG. 2: (a) Correlation length £ (on a logarithmic scale) plot- 
ted vs. T for Ising strips of width D = 10. Monte Carlo 
results (shown with error bars) were extracted for a system 
of L = 10000, recording the wavevector-dependent suscepti- 
bility X (k) = N{\M(k)\ 2 ) where M(k) = £ Sj exp(ik ■ fj) 
for k oriented in the long direction and k — k m i n = 2ir/L, 
using the formula quoted in the text. Transfer matrix re- 
sults were computed from the exact formula {Eq.(4.39)of [22j) 
for the flxoo system. Broken curve shows the approxima- 
tion £ » exp(Da /T) where a is the exactly known interfa- 
cial tension of the Ising model. The value of £ at T c [211 ]. 
£ c = (4D/7r), is shown as a dot. (b) Weight IV of the central 
peak for strips of width D = 10 plotted vs. temperature for 
L = 1000, 480, 240, 120, 80, 60 from the left to the right. The 
symbols indicate £ = L, L/3 and L/9, respectively. Insert 
shows a plot of To(L,D) vs. L (note logarithmic scale) for 
two choices of D. 



nucleation of interfaces is sufficiently easy. 

In order to show that these results carry over to flu- 
ids confined in long cylindrical pores, we have studied 
the Asakura-Oosawa (AO) model [37J of colloid-polymer- 
mixtures. The latter system is attractive for experiments: 
the large colloid size renders effects of the atomistic corru- 
gation of pore walls negligible, and facilitates observation 
of wetting layers and interfaces |38( . We describe colloids 
as hard spheres of radius r c = 1, and polymers as soft 
spheres of radius r p = 0.8. Polymer-colloid overlap (as 
well as colloid-colloid overlap) is strictly forbidden, while 
polymers can overlap with no energy cost. The phase 
diagram of this model in the bulk and for thin films has 
already been carefully studied 39|, |40( . At the cylinder 
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FIG. 3: (a) Barrier AF/T against nucleation of interfaces in 
Ising strips plotted vs. T. Several choices of L and D are 
shown, as indicated, (b) Barrier AF/T against nucleation 
of interfaces in the AO model confined to cylindrical pores 
of diameter D — 12 plotted vs. inverse polymer reservoir 
packing fraction 1/rjL. 



radius R = D/2 we apply a hard wall, which may over- 
lap with neither colloids nor polymers. This boundary 
condition at the surface leads to an entropic attraction 
of the colloidal particles to the wall [40j , causing the for- 
mation of a precursor of a wetting layer (a true wetting 
layer can only form in the limit D — > oo, of course [40j). 
For this model, the polymer fugacity exp(/ip/fcsT) or 
the related "polymer reservoir packing fraction" 77^ = 
(47rrp/3) exp^/ft^T) plays the role of an inverse tem- 
perature like variable, while the colloid packing fraction 
r/ c = (4irrl/3)N c /V (N c is the number of colloids in the 
system of volume V = irR L) is the order parameter 
density. Fig. QJa), as a counterpart of Fig. [TJb), shows 
P{rj c ) for various values of rf v . (The same Grand Canon- 
ical Monte Carlo methods as in [39[ are used.) One can 
clearly distinguish the crossover from an (asymmetric) 
double-peak distribution to a structure with three peaks, 
and finally a single peak, which only narrows when rL is 
close to r)p Crit = 0.766 [39(. The "phase diagram", where 
the coexisting polymer-rich and colloid-rich phases are 
estimated from the left-most to the right-most peak in 
Fig. BIa), is shown in Fig. 0Jb). Fig. [3fb) shows that 
the barrier against nucleation of interfaces across the 
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FIG. 4: (a) Distribution P(n c ) of the number of colloids 
in a cylinder of diameter D = 12 and length L = 180 
(all lengths are measured in units of the colloid radius) for 
rf p = 1.075,1.10,1.115,1.118,1.20 from top to bottom at 
(r/c) ~ 0.175. Above the plot we show a typical snapshot 
(cross section through the cylinder) at rf v — 1.10 containing 
multiple domains, (b) Phase diagram of the AO model in a 
cylindrical pore of diameter D — 12 and lengths L — 30, 60 
and 180, as indicated. The full curve is the bulk coexistence 
curve [33| . Note that the points shown near (r/ c ) ~ 0.16 to 
0.17 mark r]p (D, L) for three choices of L. 



pore strongly decreases with increasing L, and we have 
checked [34| that hysteresis disappears when the barrier 
is a few ksT, as for the Ising model. 

In summary, we have clarified the nature of phase co- 
existence between vapor and liquid phases of fluids (or 
fluid-fluid phase coexistence of mixtures) in long cylin- 
drical pores, depending on pore length L and pore di- 
ameter D. While at high temperatures the structure of 
the fluid is axially symmetric, phase separation in ax- 
ial direction sets in at the coexistence pressure when the 
correlation length (of the density fluctuations) £ grows 
to the order of D. Below the pore critical temperature 
£ measures the distance between domain walls, and at a 
much lower temperature (where £ « L/i) a second (again 
rounded) transition occurs (the pore then is cither in an 
axially homogeneous vapor-like or liquid- like state). The 



onset of adsorption hysteresis in the capillary is linked to 
this lower transition. A wetting transition (possible at 
a flat surface of a semi-infinite system) is also expected 
to be strongly rounded in narrow pores, and should not 
change the above conclusions. Our findings provide in- 
sight to understand experiments and simulations of fluids 
in pores, explaining the existence of a "hysteresis critical 
point" distinct from the pore critical point. A prediction 
that experiments could test is the decrease of the hys- 
teresis critical temperature with increasing pore length. 
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